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ABSTRACT 

Recently,  Goldberg  proposed  a  new  approach  lo  the  maximum  netv.ork 
flow  problem.  The  approach  yields  a  ver>'  simple  algorithm  nmning  in  0{n^) 
time  on  n-vertex  networks.  Incorporation  of  the  dynamic  tree  data  structure  of 
Sleator  and  Tarjan  yields  a  more  complicated  algorithm  with  a  running  time  of 
0{rvn  log  {n^ Im))  on  m-arc  networks.  Ahuja  and  Orlin  developed  a  variant  of 
Goldberg's  algorithm  that  uses  scaling  and  runs  in  0(,rm  +  n'  logL')  time  on  net- 
works with  integer  arc  capacities  bounded  by  U.  In  this  paper  we  explore  possi- 
ble improvements  to  the  Ahuja-Oriin  algorithm.  First  we  obtain  an  improved  run- 
ning lime  of  0(,nm  -t-  n^  logL'/log  logL')  by  using  a  nonconsiant  scaling  factor. 
Then  we  obtain  an  even  better  bound  of  0{run  -t-  n^Qogt/)'^)  by  combining  the 
Ahuja-Orlin  algorithm  with  the  wave  algorithm  of  Tarjan.  Fmally.  we  show  that 
the  use  of  dynamic  trees  in  the  laner  algorithm  reduces  the  running  time  to 

0{rvn  log( —  Gogt/)'  '  +  2)).  This  result  shows  that  the  combined  use  of  three 
tn 

different  techniques  results  in  speed  not  obtained  by  using  any  of  the  techniques 
alone. 
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1.  Introduction 

We  consider  algorithms  for  the  classical  maximum  network  flow  problem  [5,6,13,15,20]. 
We  formulate  the  problem  as  follows.  Let  G  =  (V.E)  be  a  directed  graph  with  vertex  set  V  and 
arc  set  E.  The  graph  G  is  iflow  network  if  it  has  two  distinct  distinguished  vertices,  a  source  s 
and  a  sink  t,  and  a  non-negative  real-valued  capacity  «(v,>v)  on  each  arc  (v,^)  €  E.  We  assume 
that  G  is  symmetric,  i.e.  (v.w)  e  £  iff  (w.v)  e  E.  We  denote  by  n,  m,  and  U  the  number  of  ver- 
tices, the  number  of  arcs,  and  the  maximum  arc  capacity,  respectively.  For  ease  in  stating  time 
bounds,  we  assume  m^n  and  t/  ^  4.  Bounds  containing  U  are  subject  to  the  assumption  that  all 
arc  capacities  are  integral.  All  logarithms  in  the  paper  are  base  two  unless  an  explicit  base  is 
given. 

A  flow  f  on  a  netwoiit  G  is  a  real-valued  function /on  the  arcs  satisfying  the  following  con- 
straints: 

f(v,w)  <  u(v,w)  for  all  (v.w)  e  E  (capacity  constraint),  (1) 

f(y,w)  =  -f(w,v)  for  all  (v.w)  e  E  (antisymmetry  constraint),  (2) 

S     /(v.**')  =  0  for  all  w  e  V-[s.t]  (conservation  constraint).  (3) 

(v.w)e£ 
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The  value  I /I  of  a  flow /is  the  nel  flow  into  the  sink: 


1/1=    I    /(v./). 

(v.Oe£ 


A  maximum  flov.-  is  a  flow  of  maximum  value.  The  maximum  flov,'  problem  is  that  of  finding  a 
maximum  flow  in  a  given  network. 

Remark.  We  assimie  that  all  arc  capacities  are  finite.  If  some  arc  capacities  are  infinite  but  no 
path  of  infinite-capacity  edges  from  j  to  r  exists,  then  each  infinite  capacity  can  be  replaced  by  the 
sum  of  the  finite  capacities,  without  affecting  the  problem.  D 

The  maximum  flow  problem  has  a  long  and  rich  history,  and  a  series  of  faster  and  faster 
algorithms  for  the  problem  has  been  developed.  (See  [10]  for  a  brief  survey.)  The  previously 
fastest  knowTi  algorithms  are  that  of  Goldberg  and  Tarjan  [8,10],  with  a  running  time  of 
0{,nm  log  (n^/m)),  and  that  of  Ahuja  and  Orlin  [1],  with  a  running  time  of  0{,nm  +  n}  logt/). 
Both  of  these  algorithms  are  refinements  of  a  generic  method  proposed  by  Goldberg  [8],  which 
we  shall  call  the  prefknv  algorithm.  For  networks  with  m  =  fi(n^),  the  Goldberg-Tarjan  bound  is 
O(n^),  which  matches  the  bound  of  several  eariier  algorithms  [12,14,16,21].  For  networks  with 
m  =  0(n^"^)  for  some  constant  e  >  0,  the  Goldberg-Tarjan  bound  is  0(nm  \ogn),  which  matches 
the  bound  of  the  earlier  Sleator-Tarjan  algorithm  [17,18].  Under  the  similarity  assumption  [7], 
namely  U  =  0(n*)  for  some  constant  k,  the  Ahuja-Oriin  bound  beats  the  Goldberg-Tarjan  bound 
unless  m  =  0(n)orm  =  Cl{n^). 

The  Goldberg-Tarjan  and  Ahuja-Orlin  algorithms  obtain  their  speed  from  two  different 
techniques.  The  former  uses  a  sophisticated  data  structure,  the  dynamic  tree  structure  of  Sleator 
and  Taijan  [18,19,20],  whereas  the  latter  uses  scaling.  In  this  paper  we  explore  improvements  in 
the  Ahuja-Oriin  algorithm  obtainable  by  incorporating  other  ideas,  including  the  use  of  dynamic 
trees.  We  begin  in  Section  2  by  reviewing  the  generic  preflow  algorithm  [8,9,10].  In  Section  3, 
we  develop  a  version  of  the  Ahuja-Oriin  algorithm  that  uses  a  stack-based  vertex  selection  rule 
and  a  nonconstanl  scaling  factor  to  obtain  a  time  bound  of  0(nm  +  n^  logt//log  logt/).  In  Sec- 
tion 4,  we  describe  an  even  faster  variant,  which  uses  a  coristani  scaling  factor  but  combines  the 
scaling  idea  with  the  wave  algorithm  of  Tarjan  [21].  This  algorithm  has  a  running  time  of 
0(nm  -»-  fi^  QogU)^'^).  In  Section  5  we  add  dynamic  trees  to  the  method  of  Seaion  4,  thereby 

obtaining  a  running  time  of  0(nm  log(  —  Oogt/)"^  +2).  The  results  in  Sections  3-5  rely  on  the 

m 

assumption  that  pointer  manipulations  and  arithmetic  operations  on  integers  of  magnitude  U  take 
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0(1)  time.  In  Section  6  we  consider  the  effect  on  our  time  bounds  of  a  weaker  assumption, 
namely  that  arithmetic  on  integers  of  magnitude  n  takes  0(1)  time.  We  conclude  in  Section  7 
with  some  final  remarks. 

2.  The  Preflow  Algorithm 

In  contrast  to  the  classical  augmenting  path  method  of  Ford  and  Fulkerson  [6],  which  moves 
flow  along  an  entire  path  from  ;  to  f  at  once,  the  preflow  method  moves  flow  along  a  single  arc  at 
a  time.  The  key  concept  underlying  the  algorithm  is  that  of  a  preflow,  introduced  by  Karzanov 
[12].  A  preflow /is  a  real-valued  function  on  the  arcs  satisfying  constraints  (1),  (2),  and  a  relaxa- 
tion of  (3).  For  any  vertex  w.  let  the  flow  excess  of  w  be  e(w)  =     "^     /(v.w).  The  required 

(v,w)€  E 

constraint  is  the  following: 

e(w)  ^  0  for  all  w  €  V-{s}  (nonnegativity  constraint).  (4) 

We  call  a  vertex  v  active  if  v^  /  and  e(y)  >  0.  Observe  thai  the  nonnegativity  constraint  implies 
thatc(5)<0. 

The  residual  capacity  of  an  arc  (v,>v)  with  respect  to  a  preflow  /  is 
uj{v,w)  =  u(v,w)-f(v,w).  An  arc  is  saturated  if  ujiv,w)  =  0  and  unsaturated  otherwise.  (The 
capacity  constraint  implies  that  any  unsaturated  arc  (v.w)  has  u/v.w)  >  0.) 

The  preflow  algorithm  maintains  a  preflow  and  moves  flow  from  aaive  vertices  through 
unsaturated  arcs  toward  the  sink,  along  paths  estimated  to  contain  as  few  arcs  as  possible.  Excess 
flow  that  cannot  be  moved  to  the  sink  is  returned  to  the  source,  also  along  estimated  shortest 
paths.  Eventually  the  preflow  becomes  a  flow,  which  is  a  maximum  flow. 

As  an  estimate  of  path  lengths,  the  algorithm  uses  a  valid  labeling,  which  is  a  function  d 
from  the  vertices  to  the  nonnegative  integers  such  that  dis)  =  n,  d{t)  =  0,  and  d(v)  <  d(w)  +  1  for 
every  unsaturated  arc  (v.h').  A  proof  by  induction  shows  that,  for  any  valid  labeling  d, 
d(v)  ^  min  [df{v,s)  +  n,dj{v,t)],  where  dj{v,w)  is  the  minimum  number  of  arcs  on  a  path  from  v 
to  w  consisting  of  arcs  unsaturated  with  respect  to  the  flow  /.  We  call  an  arc  (v,>v)  eligible  if 
(v,w)  is  unsaturated  and  d{v)  =  d(w)  +  1. 


.4- 


The  algorilhm  begins  with  an  initial  preflow/and  valid  labeling  J  defined  as  follows: 


/(V.H')=' 


u{v,w)  if  V  =s, 
-u(w,v)  if  w  =  s, 
0  if  V  *  J  and  w  ^s. 


d(v)  =  m\n  [df(v,s)  +  n,djiv,t)]. 
The  algorithm  consists  of  repeating  the  following  two  steps,  in  any  order,  until  no  vertex  is  ac- 


tive: 


Push  (y,w). 

Applicability:  Vertex  v  is  active  and  arc  (v,w)  is  eligible. 

Action:  Increase /(v.w)  by  min  [e(v),uj{v,w)].  The  push  is  saturating  if  (v.w)  is 
saturated  after  the  push  and  nonsaturating  otherwise. 

Relabel  (v). 

Applicability:  Vertex  v  is  active  and  no  arc  (y,w)  is  eligible. 
Action:  Replace  d(v)  by  min  {d(w)  +1  I  (v.w)  is  unsaturated). 

When  the  algorithm  terminates, /is  a  maximum  flow.  Goldberg  arxl  Tarjan  derived  the  fol- 
lowing bounds  on  the  number  of  steps  required  by  the  algorithm: 

Lemma  2.1  [10].  Relabeling  a  vertex  v  strictly  increases  d{v).  No  vertex  label  exceeds  2n-l. 
The  total  number  of  relabelings  is  0(,n^). 

Lemma  2.2  [10].  There  are  at  most  0{nm)  saturating  pushes  and  at  most  0{n^m)  nonsaturating 
pushes. 

Efficient  implementations  of  the  above  algorithm  require  a  mechanism  for  selecting  pushing 
and  relabeling  steps  to  perform.  Goldberg  arxl  Tarjan  proposed  the  following  method.  For  each 
vertex,  construct  a  (fixed)  list  A{y)  of  the  arcs  out  of  v.  Designate  one  of  these  arcs,  initially  the 
first  on  the  Lst,  as  the  current  arc  out  of  v.  To  execute  the  algorithm,  repeat  the  following  step 
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until  there  no  active  vertices: 

PushJRelabel  (v). 

Applicability:  Vertex  v  is  active. 

Action:  If  the  current  arc  (v.h-)  of  v  is  eligible,  perform  push(y,w).  Othenvise.  if  (v.w) 

is  not  the  last  arc  on  i4(v),  make  the  next  arc  after  (v.w)  the  current  one.  Other- 
wise, perform  relabel  (v)  and  make  the  first  arc  on  /4(v)  the  current  one. 

With  this  implementation,  the  algorithm  runs  in  0(jvn)  time  plus  0{\)  time  per  nonsaturat- 
ing  push.  This  gives  an  0{,n}m)  time  bound  for  any  order  of  selecting  vertices  for  push/relabel 
steps.  Making  the  algorithm  faster  requires  reducing  the  time  spent  on  rwnsaturating  pushes.  The 
number  of  such  pushes  can  be  reduced  by  selecting  vertices  for  push/relabel  steps  carefully. 
Goldberg  and  Tarjan  showed  that  FIFO  selection  (first  active,  first  selected)  reduces  the  number 
of  rxjnsaturating  pushes  to  0{n}).  Cheriyan  and  Maheshwari  [3]  showed  that  highest  label  selec- 
tion (always  pushing  flow  from  a  vertex  uiih  highest  label)  reduces  the  number  of  nonsaturating 
pushes  to  0{n}m^'^).  Ahuja  and  Orlin  proposed  a  third  selection  rule,  which  we  discuss  in  the 
next  section. 

3.  The  Scaling  Algotithm 

The  intuitive  idea  behind  the  Ahuja-Oriin  algorithm,  henceforth  called  the  scaling  algo- 
rithm, is  to  move  large  amounts  of  flow  when  possible.  The  same  idea  is  behind  the  maximum 
capacity  augmenting  path  method  of  Edmonds  and  Karp  [4]  and  the  capacity  scaling  algorithm  of 
Gabow  [7].  One  way  to  apply  this  idea  to  the  preflow  algorithm  is  to  always  push  flow  from  a 
vertex  of  large  excess  to  a  vertex  of  small  excess,  or  to  the  sink.  The  effect  of  this  is  to  reduce  the 
maximum  excess  at  a  rapid  rate. 

Making  this  method  precise  requires  specifying  when  an  excess  is  large  and  when  small. 
For  this  purpose  the  scaling  algorithm  uses  an  excess  bound  A  arxl  an  integer  scaling  factor  kt2. 
A  vertex  v  is  said  to  have  large  excess  if  its  excess  exceeds  A/k  and  small  excess  otherwise.  As 
the  algorithm  proceeds,  k  remains  fixed,  but  A  periodically  decreases.  Initially,  A  is  the  smallest 
power  of  k  such  that  A  S  I/.  The  algorithm  maintains  the  invariant  that  e{v)  S  A  for  every  active 
vertex  v.  This  requires  changing  the  pushing  step  to  the  following. 

Push  (v.w). 

Applicability:  Vertex  v  is  active  and  arc  (v.w)  is  eligible. 

Action:  If  w  *r,  increase /(v.w)  by  min  (f(v).«/(v,H'),  A-*(h')).  Otherwise  (w  =  r), 
increase /(v.w)  by  min  [e{v),Uf{v,w)). 


The  algorithm  consists  of  a  number  of  scaling  phases,  during  each  of  which  A  remains  con- 
stant. A  phase  consists  of  repeating  push/relabel  steps,  using  the  following  selection  rule,  until 
no  active  vertex  has  large  excess,  and  then  replacing  A  by  A/it.  The  algorithm  terminates  when 
there  are  no  active  vertices. 

Large  excess,  smallest  label  selection:  Apply  a  push/relabel  step  to  an  active  vertex  v  of  large 
excess;  among  such  vertices,  choose  one  of  smallest  label. 

If  the  edge  capacities  are  integers,  the  algorithm  terminates  after  at  most  [  log*  U  +  \} 
phases:   after  [  log*  t/+  ij   phases,  A  <  1,  which  implies  that  /is  a  flow,  since  the  algorithm 
mainuins  integrality  of  vertex  excesses.  Ahuja  and  Orlin  derived  a  bound  of  0(kn^  log*  U)  on 
the  total  number  of  nonsaturating  pushes.  We  repeat  the  analysis  here,  since  it  provides  motiva- 
tion for  our  first  modification  of  the  algorithm. 

Lemma  3.1   [1].    The  total  number  of  nonsaturating  pushes  in  the  scaling  algorithm   is 
OC/k/j^OogilZ  +  l)). 

Proof.   Consider  the  function  <!)=    J^  e(v)rf(v)/A.   We  call  O  the  potential  of  the  current 

V  acbve 

preflow  /  and  labeling  d.  Since  0  <  e(v)/A  <  1  and  0  <  d(v)  <  2n  for  every  active  vertex  v, 
0  <  <1>  <  2«^  throughout  the  algorithm.  Every  pushing  step  decreases  <I>.  A  nonsaturating  push- 
ing step  decreases  <1>  by  at  least  l/k.  since  the  push  is  from  a  vertex  v  with  excess  more  than  A/k 
to  a  vertex  w  with  d(w)  =  d(v)  -1,  and  e(w)  ^  A/k  or  w  =  r.  The  value  of  <I>  can  increase  only 
during  a  relabeling  or  when  A  changes.  A  relabeling  of  a  vertex  v  increases  <1>  by  at  most  the 
amount  ^(v)  increases.  Thus  the  total  increase  in  A  due  to  relabelings,  over  the  entire  algorithm, 
is  at  most  2n^.  When  A  changes.  <I>  increases  by  a  factor  of )!:,  to  at  most  2n^.  This  happens  at 
most  [  logi  U  +  \j  times.  Thus  the  total  increase  in  <I>  over  the  entire  algorithm  is  at  most 
2n^  I  log*  U  +  2}  .  The  total  number  of  nonsaturating  pushes  is  at  most  k  times  the  sum  of  the 
initial  value  of  <I>  and  the  total  increase  in  C).  This  is  at  most  2kn^  [  log^  t/  +  3j  .  D 

Choosing  /t  to  be  a  constant  independent  of  n  gives  a  total  time  bound  of  0(nm  -«-  n^  log  U) 
for  the  scaling  algorithm,  given  an  efficient  implementation  of  the  vertex  selection  rule.  One  way 
to  implement  the  rule  is  to  maintain  an  array  of  sets  indexed  by  vertex  label,  each  set  containing 
all  large  excess  vertices  with  the  corresponding  label,  and  to  maintain  a  pointer  to  the  nonempty 
set  of  smallest  index.  The  total  time  needed  to  maintain  this  structure  is  0(nm  +  n^  log  U). 

Remark.  The  bound  for  the  scaling  algorithm  can  be  improved  slightly  if  it  is  measured  in  terms 
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of  a  different  parameter.    Let  (/*=4+     J^     u(s,v)/n.    Then  the  bound  on  nonsaturating 

pushes  can  be  reduced  to  O  (n^  log  U')  and  the  overall  time  bound  to  O  (nm  +  n}  log  (/*).  This 
improvement  is  analogous  to  the  bound  Edmonds  and  Karp  derived  for  their  capacity-scaling 
transportation  algorithm  [4].  The  argument  is  as  follows.  The  algorithm  maintains  the  invariant 
that  the  total  excess  on  active  vertices  is  at  most  nU* .  Let  phase  j  be  the  first  phase  such  that 
A  ^  {/*.  Then  the  total  iixrcase  in  <I>  due  to  phase  changes  up  to  and  including  the  change  from 

phase  y  -  1  to  phase  ;  is  at  most  n  2)  2n/2'~*  =  0(n  ).  The  total  increase  in  <I>  due  to  later  phase 
changes  is  O  (n^  logt/*).  D 

Having  described  the  scaling  algorithm,  we  consider  the  question  of  whether  its  running 
time  can  be  improved  by  reducing  the  number  of  nonsaturating  pushes.  The  proof  of  Lemma  3.1 
bounds  the  number  of  nonsaturating  pushes  by  estimating  the  total  increase  in  the  potential  <I>. 
Observe  that  there  is  an  imbalance  in  this  estimate:  0{n^  log*  U)  of  the  increase  is  due  to  phase 
changes,  whereas  only  O(n^)  is  due  to  relabelings.  Our  plan  is  to  improve  this  estimate  by 
decreasing  the  contribution  of  the  phase  changes,  at  the  cost  of  increasing  the  contribution  of  the 
relabelings.  Making  this  plan  work  requires  changing  the  algorithm. 

We  use  a  nonconstant  scale  factor  and  a  slightly  more  elaborate  method  of  vertex  selection. 
We  make  use  of  the  stack-push/rtldbcl  step  defined  below,  which  performs  a  sequence  of  push 
and  relabel  steps  using  a  stack.  The  stack  provides  an  alternative  way  of  avoiding  pushes  to 
large-€xcess  vertices  (other  than  r). 

Stack-Push/ReUibel(r). 

Applicability:  Vertex  r  is  active. 

Action:  Initialize  a  stack  5  to  contain  r  Repeat  the  following  step  until  5  is  empty: 

Stack  Step.  Let  v  be  the  top  vertex  on  S  and  let  (v,^)  be  the  current  arc  out  of  v. 
Apply  the  appropriate  one  of  the  following  cases: 

Case  1:  (v,^)  is  not  eligible. 

Case  la:  {v,-w)  is  not  last  on  A(v).  Replace  (v,**-)  as  the  current  arc  out  of  v 

by  the  next  arc  on  y4(v). 
Case  lb:  {y,w)  is  last  onA(v).  Relabel  v  and  pop  it  from  5.  Replace  (v.w)  as 

the  current  arc  out  of  v  by  the  first  arc  on  A(v). 
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Case  2:  (v.h')  is  eligible. 

Case  2a:  ^Cw)  >  A/2  and  w*^t.  Push  k-  onto  S. 

Case  2b:  e{w)  <  till  or  w  =  r.  Perform  push{\,-w)  (modified  as  at  the  begin- 
ning of  this  section  to  maintain  e(yv)  <  A  if  w  *  r).  If  e{y)  =  0,  pop  v 
from  S. 

Some  remarks  about  stack-push! relabel  are  in  order.  Let  us  call  a  nonsaturating  push  big  if 
it  moves  at  least  A/2  units  of  flow  and  little  otherwise.  During  an  execution  of  stack- 
push/relabel,  every  vertex  pushed  onto  S,  except  possibly  the  first,  has  an  excess  of  at  least  A/2 
when  it  is  added  to  S.  A  vertex  v  can  only  be  popped  from  S  after  it  is  relabeled  or  its  excess  is 
reduced  to  zero.  Of  the  pushes  from  v  while  v  is  on  S,  at  most  two  are  nonsaturating,  only  the  last 
of  which  can  be  little. 

Our  variant  of  the  scaling  algorithm,  called  the  stack  scaling  algorithm,  consists  of  phases 
just  as  in  the  scaling  algorithm.  A  phase  consists  of  repeatedly  applying  the  stack-pushi relabel 
step  to  a  large-excess  active  vertex  of  highest  label;  a  phase  ends  when  there  are  no  large-excess 
active  vertices. 

Lemma  3.2.  The  total  number  of  nonsaturating  pushes  made  by  the  stack  scaling  algorithm  is 
Oikn'^-t-n^QogtU  +  l)). 

Proof.  To  bound  the  number  of  nonsaturating  pushes,  we  use  an  argument  like  the  proof  of 
Lemma  3.1,  but  with  two  potentials  instead  of  one.  The  first  potential  is  that  of  Lemma  3.1, 
namely  0=    J^    e(,v)d(v)/A.  By  the  analysis  in  the  proof  of  Lemma  3.1,  every  push  decreases 

V  tctive 

<I),  the  total  increase  in  <1>  over  all  phases  is  0(n^  logi  U),  and  the  difference  between  the  initial 
and  the  final  values  of  C)  is  O(n^).  Each  big  push  moves  at  least  A/2  imits  of  flow,  and  hence 
decreases  <I)  by  at  least  1/2.  Thus  the  number  of  big  pushes  is  O(n^(\ogt  U  +  1)). 

To  count  little  pushes,  we  divide  them  into  two  kinds,  those  that  result  in  an  empty  stack  5, 
called  emptying  pushes,  and  those  that  do  not,  called  nonemptying  pushes.  A  nonemptying  push 
from  a  vertex  v  is  such  that  e(v)  was  at  least  A/2  when  v  was  added  to  5,  and  the  push  results  in 
e(v)  decreasing  to  zero.  We  can  charge  such  a  push  against  the  cimiulative  decrease  of  at  least 
1/2  in  <1>  resulting  from  moving  the  original  excess  on  v  to  vertices  of  smaller  label.  Hence  there 
can  be  only  0{n^ Qogt  (/+!))  rwnemptying  pushes. 

An  emptying  push  from  a  vertex  v  can  be  associated  with  a  decrease  of  at  least  l/k  in  <t>, 
namely  the  drop  in  <t)  caused  by  the  movement  of  the  original  excess  on  v,  which  is  at  least  A/k, 


to  smaller  labeled  vertices.  But  using  this  drop  gives  a  bound  on  the  number  of  emptying  pushes 
of  only  Oikn^Qogt  U  +  1)).  We  count  emptying  pushes  more  carefully  by  using  a  second  fwten- 
tial,  <1>2.  The  definition  of  <1>2  involves  two  parameters,  an  integer  /  and  a  set  P.  The  value  of  /  is 
equal  to  the  minimum  of  2n  and  the  smallest  label  of  a  vertex  added  to  S  while  S  was  empty  dur- 
ing the  current  phase.  Observe  that  /  =  2n  at  the  beginning  of  a  phase  and  that  /  is  nonincreasing 
during  a  phase.  The  set  P  consists  of  all  vertices  with  label  greater  than  /  and  all  vertices  with 
label  equal  to  /  from  which  an  emptying  push  has  been  made  during  the  current  phase.  Observe 
that  P  is  empty  at  the  beginning  of  a  phase  and  P  never  loses  a  vertex  during  a  phase.  The 
definition  of  <I>2  is 


<I>2=        Z        e(v)(rf(vH+l)/A.  (lf/»  =  0.  <D2  =  0.) 

v€  P:e(v)>0 


Observe  that  0  <  O2  ^  2n^.  Any  emptying  push  can  be  associated  either  with  an  addition  of 
a  vertex  to  /*  or  with  a  decrease  in  <1)2  of  at  least  1  /k.  (If  the  push  is  from  v,  either  v  4.  P  when  v 
is  added  to  S  but  v  e  P  after  the  push,  orv  e  P  when  v  is  added  to  S  and  <1>2  drops  by  at  least  \/k 
because  of  pushes  from  v  while  v  is  on  5.)  The  number  of  vertices  added  to  P  is  at  most 
n  I  logt  (/  +  Ij  over  all  phases,  and  hence  so  is  the  niunber  of  emptying  pushes  not  associated 
with  decreases  in  <I>2. 

To  bound  the  number  of  emptying  pushes  associated  with  decreases  in  <I)2,  we  bound  the 
total  increase  in  <1>2.  Increases  in  <l>2  are  due  to  relabelings  and  to  decreases  in  /.  (A  vertex  added 
to  P  because  of  an  emptying  push  has  zero  excess  and  hence  adds  nothing  to  <1>2.)  A  relabeling  of 
a  vertex  v  increases  <&2  by  at  most  the  increase  in  d(v)  plus  one;  the  "plus  one"  accounts  for  the 
fact  that  the  relabeling  may  add  v  to  F.  Thus  relabelings  contribute  at  most  4n^  to  the  growth  of 

There  are  at  most  2n  decreases  in  /  per  phase.  A  decrease  in  /  by  one  adds  at  most  n/k  to 
<I>2.  since  when  the  decrease  occurs  every  vertex  in  P  has  small  excess.  Thus  the  total  increase 
in  <^2  due  to  decreases  in  /  is  at  most  2n^  [  log*  U  +\\  Ik  over  all  [rfiases. 

The  total  nimiber  of  emptying  pushes  associated  with  decreases  in  <I>2  is  at  most  k  times  the 
total  increase  in  <I>2.  since  <I>2  is  initially  zero.  Thus  the  total  number  of  such  pushes  is  at  most 
4/bi2+2/jn  logifZ  +  lJ  .n 

As  in  the  Goldberg-Tarjan  and  Ahuja-Oriin  algorithms,  the  time  to  perform  saturating 
pushes  and  relabeling  operations  is  0{nm).  The  only  significant  remaining  issue  is  how  to 
choose  vertices  to  add  to  S  when  S  is  empty.   For  this  purpose,  we  maintain  a  data  structure 
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consisting  of  a  coUection  of  doubly  linked  lists /«f(/)  =  (' e  N :  e(i)>  A/ k  and  d(i)  =  j]  for  each 

j  €  (1,2 2n-l ).  We  also  maintain  a  pointer  to  indicate  the  largest  index  j  for  which  listij)  is 

nonempty.  Maintaining  this  structure  requires  0(1)  time  per  push  operation  plus  time  to  main- 
tain the  pointer.  The  pointer  can  only  increase  due  to  a  relabeling  (by  at  most  the  amount  of 
change  in  label)  or  due  to  a  phase  change  (by  at  most  2n).  Consequently,  the  number  of  times  the 
pointer  needs  to  be  incremented  or  decremented  is  0(n^  +  n  Oog*  U  +  1)).  The  overall  running 
lime  of  the  algorithm  is  thus  0(nm  +  kn^  +  n^  (log*  U  +  1)).  Choosing  k  ={  log  U /log  log  U] 
gives  the  following  result: 

Theorem  3.3.  The  stack  scaling  algorithm,  with  an  appropriate  choice  of  k,  nms  in 
O  {run  +  n^  log  t//loglog  U)  time. 

Remark.  This  bound  can  be  improved  to  O  (nm +  n^  log  f/' /log  log  t/*).  where 
t/*  =4  +      ^      u(s,v)/n.  The  algorithm  must  be  changed  slightly.  We  use  a  scaling  factor  of 

v:(i,v)€  E 

k  =  2  until  A  <  [/'  and  then  switch  to  it  =[  log  U*  /log  log  U"\  for  the  remainder  of  the  algo- 
rithm. When  the  switch  occurs,  we  increase  A  (if  necessary)  to  be  a  power  of  it.  The  analysis  of 
the  initial  phases  (those  in  which  A  >  t/')  is  the  same  as  that  in  the  remark  following  Lemma  3.1. 
These  phases  account  for  only  O(n^)  nonsaturating  pushes.  The  analysis  of  the  remaining  phases 
is  as  above  with  U  replaced  by  U* .O 


4.  The  Wave  Scaling  Algorithm 

Another  way  to  reduce  the  number  of  nonsaturating  pushes  in  the  scaling  algorithm  is  to 
keep  track  of  the  total  excess,  defined  to  be  the  sum  of  the  excesses  of  all  active  vertices.  The  key 
observation  is  that  if  the  total  excess  is  sufficiently  large,  the  algorithm  can  make  significant  pro- 
gress by  applying  slack-pushJ relabel  to  each  active  vertex  in  turn,  processing  vertices  in  topologi- 
cal order  with  respea  to  the  set  of  eligible  arcs.  Even  though  some  vertices  of  very  small  excess 
are  processed,  the  overall  benefits  of  this  approach  yield  an  even  better  running  time  than  that  of 
the  Slack  scaling  algorithm.  The  idea  of  processing  vertices  in  topological  order  originated  in  the 
wave  algorithm  of  Taijan  [20,21];  therefore  we  call  the  new  algorithm  the  wave  scaling  algo- 
rithm. 

The  wave  scaling  algorithm  seems  to  derive  no  benefit  from  using  a  nonconstant  scaling 
factor,  therefore  we  fix  it  =  2.  The  algorithm  uses  another  parameter  /  S  1  whose  exact  value  we 
shall  choose  later.  A  phase  of  the  wave  scaling  algorithm  consists  of  two  pans.  First,  the  wave 
step  below  is  repeated  until  the  total  excess  is  less  than  n  A//.  Then  stack-pusk' relabel  steps  are 
applied  to  large-excess  active  vertices  in  any  order  until  there  are  no  large-excess  active  vertices. 
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Wave;  Construct  a  list  L  of  the  vertices  in  V  -  [s,t]  in  topological  order  with  respect  to  the  set  of 
eligible  arcs.  (Ordering  L  in  nonincreasing  order  by  vertex  label  suffices.)  For  each  ver- 
tex V  on  L,  if  V  is  active  perfonn  stack-pushlrelabel{v). 

Observe  that  the  total  excess  is  a  nonincreasing  function  of  time.  Constructing  L  at  the 
beginning  of  a  wave  takes  0{n)  time  using  a  radix  sort  by  vertex  label.  The  time  spent  during  a 
single  wave  is  0(n)  plus  time  for  the  relabelings  and  pushes.  Thus  the  total  time  required  by  the 
wave  scaling  algorithm  is  0(nm)  plus  0(n)  per  wave  plus  0(1)  per  nonsaturating  push.  We 
complete  the  running  time  analysis  with  two  lemmas. 

Lemma  4.1.  The  number  of  nonsaturating  pushes  done  by  the  wave  scaling  algorithm  is 
0(n^  -K  (n^//)logt/)  plus  0(n)  per  wave. 

Proof.  We  coimt  big  pushes  (those  that  move  at  least  A/2  units  of  flow)  and  little  pushes 
separately.  The  analysis  in  the  proof  of  Lemma  3. 1  applies  to  bound  the  number  of  big  pushes, 
with  the  following  improvement.  At  the  end  of  a  phase,  the  total  excess  is  less  than  nAll.  Thus 
the  potential  at  the  beginning  of  the  next  phase  is  at  most  4n^/l.  (The  new  value  of  A  is  half  of 
the  old;  each  vertex  label  is  at  most  2n.)  Thus  the  sum  of  the  increases  in  <1>  caused  by  changes  in 
A  is  Oi(,n^/l)logU).  which  implies  that  the  number  of  big  pushes  is  0(n^  +  (/i^//)log  U).  The 
same  argument  gives  the  same  bound  on  the  number  of  little  pushes  from  vertices  that  have  large 
excess  when  added  to  the  stack  5,  since  as  in  the  pnx)f  of  Lemma  3.2  each  such  push  can  be 
charged  against  a  drop  of  at  least  1/2  in  <I>.  The  remaining  litUe  pushes  consist  of  at  most  one  per 
vertex  to  which  stack-push! relabel  is  applied  during  waves,  totaling  at  most  n  per  wave.  D 

Lemma  4 2.  The  number  of  waves  in  the  wave  scaling  algorithm  is  O  (min(n/  +  logf/,n^)). 

Proof.  Consider  any  wave  except  the  last  in  a  phase.  At  the  end  of  the  wave,  the  total  excess  is 
at  least  nA//.  The  only  way  for  excess  to  remain  on  a  vertex  x  at  the  end  of  a  wave  is  for  x  to 
have  been  relabeled  during  the  wave;  once  x  is  relabeled,  e{,x)  remains  constant  until  the  end  of 
the  wave.  The  maximum  excess  on  any  single  vertex  is  A.  It  follows  that  at  least  nil  relabelings 
must  have  occurred  during  the  wave.  Since  the  total  number  of  relabelings  is  0{n^),  the  total 
number  of  waves  is  0{,nl  +  \ogU).  Furthermore,  a  wave  during  which  no  relabeling  occurs 
causes  all  vertices  to  become  inactive,  and  hence  terminates  the  entire  algorithm.  Thus  the 
numberof  waves  is  O(n^).  D 

Theorem  43.  The  mnning  time  of  the  wave  scaling  algorithm  is  0(,nm  +  n^Oogt/)"^)  if  /  is 
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chosen  equal  lo  Qogt/)"^. 

Proof.  By  Lemmas  4.1  and  4.2,  the  running  time  of  the  wave  scaling  algorithm  is 
0(n/n -t-(n^//)logt/ -^  min  (n^/ -^  nlogt/.n^)).  Choosing  l  =  (\ogU)^'^  gives  a  bound  of 
OCwn -t- n^Oogt/)"^  +  min{«logt/.«^)).  But  min{/jlog(;.n^)  ^n^Gogt/)"^.  since 
n\ogU  S  n^Oogt/)"^  implies  Gogt/)"^  ^  n.  D 

Remark.  The  bound  in  Theorem  4.3  can  be  reduced  lo  C>(n/n -i- /i^Oogt/*)"^).  where 
t/*=4-»-     5^    M(j,v)/n,  by  doing  waves  only  after  A  <  t/*.  The  analysis  of  the  early  phases 

(j.v)€  £ 

(A  >  U*)  is  the  same  as  that  in  the  remark  following  Lemma  3.L  The  analysis  of  the  later  phases 
(A  <,  V')  is  as  above  with  U  replaced  by  U* .  D 


5.  Use  of  Dynamic  Trees 

The  approach  taken  in  Sections  3  and  4  was  to  reduce  the  total  number  of  pushes.  An 
Orthogonal  approach  is  to  reduce  the  total  time  of  the  pushes  without  necessarily  reducing  their 
number.  This  can  be  done  by  using  the  dynamic  tree  data  structure  of  Sleator  and  Tarjan 
[18,19,20].  We  conjecture  that,  given  a  version  of  the  preflow  algorithm  with  a  bound  ofp^  run 

on  the  total  number  of  pushes,  the  running  time  can  be  reduced  from  0{p)  to  0{nm  log(-^  -i- 1)) 

run 

by  using  dynamic  trees.  Although  we  do  not  know  how  to  prove  a  general  theorem  to  this  effect, 
we  have  been  able  to  obtain  such  a  result  for  each  version  of  the  preflow  algorithm  that  we  have 
considered.  As  an  example,  the  0(,nm  log  (n^/m))  time  bound  of  Goldberg  and  Tarjan  results 
from  using  dynamic  trees  with  the  FIFO  selection  rule;  the  bound  on  the  number  of  pushes  in  this 
case  is  Oin^).  In  this  section  we  shall  show  that  the  same  idea  applies  to  the  wave  scaling  algo- 
rithm of  Section  4,  resulting  in  a  time  bound  of  Oinm  log( —  (log  l/)"^  +  2)).  This  idea  can  also 

m 

be     applied     to     the     stack     scaling     algorithm     of    Section     3,     giving     a    bound     of 

O  (nm  log  ( — ; — f — —  -f  2));  we  omit  a  description  of  the  latter  result  since  the  former  bound  is 
m  log  logt/ 

better. 

The  dynamic  tree  data  structure  allows  the  maintenance  of  a  collection  of  vertex-disjoint 
rooted  trees,  each  arc  of  which  has  an  associated  real  value.  We  regard  each  tree  arc  as  being 
directed  from  child  to  parent,  and  we  regard  every  vertex  as  being  both  an  ancestor  and  a  descen- 
dant of  itself.  The  data  structure  supports  the  following  seven  operations: 
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find-root  (v):       Find  and  return  the  root  of  the  tree  containing  vertex  v. 

find-size  (v):        Find  and  return  the  number  of  vertices  in  the  tree  containing  vertex  v. 

find-value  (v):  Find  and  return  the  value  of  the  tree  arc  leaving  v.  If  v  is  a  tree  root,  the  value 
returned  is  infinity. 

find-min  (v):  Find  and  return  the  ancestor  w  of  v  with  minimum  find-valueiyv).  In  case  of  a 
tie.  choose  the  vertex  w  closest  to  the  tree  root. 

change-value  (v.x):       Add  real  number  j:  to  the  value  of  every  arc  along  the  path  from  v  to  the 
root  of  its  tree. 

link  (v.H'.x):  Combine  the  trees  containing  v  and  w  by  making  w  the  parent  of  v  and  giving 
the  arc  (v.w)  the  value  x.  This  operation  does  nothing  if  v  and  w  are  in  the  same 
tree  or  if  v  is  not  a  tree  root. 

cut  (v):  Break  the  tree  containing  v  into  two  trees  by  deleting  the  arc  joining  v  to  its 

parent;  return  the  value  of  the  deleted  arc.  This  operation  breaks  no  arc  and 
returns  infinity  if  v  is  a  tree  root. 

A  sequence  of  h  tree  operations,  starting  with  an  initial  collection  of  singleton  trees,  takes 
0{h\og  (z  -(■  1))  time  if  z  is  the  maximum  tree  size  [10.18,1920]. 

In  the  network  flow  application,  the  dynamic  tree  arcs  are  a  subset  of  the  current  arcs  out  of 
the  vertices.  Only  eligible  arcs  are  tree  arcs.  The  value  of  a  tree  arc  is  its  residual  capacity.  The 
dynamic  tree  data  structure  allows  flow  to  be  moved  along  an  entire  path  at  once,  rather  than 
along  a  single  arc  at  a  time.  We  shall  describe  a  version  of  the  wave  scaling  algorithm,  which  we 
call  the  tree  scaling  algorithm,  that  uses  this  capability.  Two  parameters  govern  the  behavior  of 
the  algorithm,  a  variable  bound  A  on  the  maximum  excess  at  an  active  vertex,  and  a  fixed  bound 
z,  1  S  z  ^  n.  on  the  maximum  size  of  a  dynamic  tree.  The  algorithm  is  identical  to  the  wave  scal- 
ing algorithm  except  that  it  uses  the  following  step  in  place  ol  stack-pushJ relabel: 

Tree-PushJRelabel{r) 

Applicability:  Vertex  r  is  active. 

Action:  Initialize  a  stack  S  to  contain  r.  Repeat  the  following  step  until  S  is  empty. 
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Stack  Step.  Let  v  be  the  top  vertex  on  S  and  let  (v,h)  be  the  current  arc  out  of  S.  Apply 
the  appropriate  one  of  the  follov.-ing  cases: 

Case  1 :  (v.h)  is  not  eligible. 

Case  la:  (v,h)  is  not  last  on  A  (v).  Replace  (v,k)  as  the  current  arc  out  of  v 

by  the  next  arc  on  y4(v). 
Case  lb:  (y\V')islastonA{y).  Relabel  v  and  pop  it  from  5.  Replace  (v,h) 

as  the  current  arc  out  of  v  by  the  first  arc  on  A  (v).  For  every  tree  arc 

O'.v),  perform  cutiy). 

Case  2:  {v,w)  is  eligible.  Ltlx=  find- root(w). 

Case  2a:  e{x)  >  t^ll  andx  *■  t.  Push  x  onto  S. 

Case2b:  e(x)<Al2orx  =  t.Ltlt=  min{e{v),Uf{v,w),find-min(\^,)}.  Lei 

6  =  eifx  =  r,  6  =  min{E,A-e(x))  if  x  *t.  Send  6  uruts  of  flow  from  v 
to  X  by  increasing  /(v,h)  by  6  and  performing  change-value  (h,-5). 
(This  is  called  a  tree  push  from  v  to  x.  The  tree  push  is 
scuurating  if  S  =  min  [uji\\w),find-min  (w)]  bcfoTZ  the  push  and 
nonsaturating  otherwise.)  While  find-value  (find-min  (w))  =  0, 
perform  cut{find-min{w)).  If  e{y)  =  0,  pop  v  from  S,  and  if  in  addition 
u/v.H')  >  0  and  find-size{yv)  +  find-size(v)  <  z,  perform  link 
(v,w,uj(v,w)). 

The  tree  scaling  algorithm  stores  flow  in  two  different  ways:  explicitly  for  edges  that  are  not 
dynamic  tree  arcs  and  implicitly  (in  the  dynamic  tree  data  structure)  for  arcs  that  are  dynamic  tree 
arcs.  After  each  cut,  the  flow  on  the  arc  cut  must  be  restored  to  its  correa  current  value.  In  addi- 
tion, when  the  algorithm  terminates,  the  correct  flow  value  on  each  remaining  tree  arc  must  be 
computed.  For  arcs  cut  during  the  computation,  the  desired  flow  values  are  returned  by  the 
corresponding  cut  operations.  Computing  correct  flows  on  termination  can  be  done  using  at  most 
n  find-value  operations.  We  have  omitted  these  bookkeeping  steps  from  our  description  of  the 
algorithm. 

The  algorithm  maintains  the  following  invariants:  every  active  vertex  is  a  tree  root,  every 
tree  arc  is  eligible,  no  excess  exceeds  A,  and  no  tree  size  exceeds  z.  Let  us  call  a  nonsaturating 
tree  push  big  if  it  moves  at  least  A/2  units  of  flow,  and  little  otherwise.  Of  the  pushes  from  a 
given  vertex  v  while  v  is  on  S,  at  most  two  are  nonsaturating  and  at  most  one  (the  last)  is  small. 
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The  tree  scaling  algorithm  is  a  variant  of  the  d>'namic  tree  algorithm  of  Goldberg  and  Tar- 
jan  [9,10].  Their  analysis  applies  to  give  the  following  result: 

Lemma  5.1  [10].  The  total  time  required  by  the  tree  scaling  algorithm  is  0{nm  log(z  +  1))  plus 
0(\og{z  +  1))  time  per  tree  push  plus  the  time  needed  to  construct  and  scan  the  list  L  during  v;ave 
steps.  The  number  of  links,  cuts,  and  saturating  tree  pushes  is  0{nm). 

Making  the  tree  scaling  algorithm  efficient  requires  careful  implementation  of  the  list  L 
used  in  wove  steps.  For  the  moment  we  shall  ignore  the  time  spent  manipulating  L.  The  remain- 
ing issue  in  analyzing  the  algorithm  is  bounding  the  number  of  nonsaturating  pushes.  To  do  this 
we  need  two  lemmas. 

Lemma  52.  The  number  of  waves  in  the  tree  scaling  algorithm  is  O  (min  [nl +  \ogU,n^)). 

Proof.  Identical  to  the  proof  of  the  same  result  for  the  wave  scaling  algorithm  (Lemma  4.2).  D 

Lemma  5.3.  The  number  of  nonsaturating  tree  pushes  done  by  the  tree  scaling  algorithm  is 
0(nm  -»-(n^//)log[/)plus  C>(rt/2)  per  wave. 

Proof.  The  potential  function  <I>  used  in  Section  3,  as  applied  in  the  proof  of  Lemma  3.3,  serves 
to  bound  the  number  of  big  tree  pushes  by  O  (n^  +  (n^//)  log  U).  We  count  little  tree  pushes  as 
follows.  Any  little  tree  push,  say  from  a  vertex  v,  reduces  ^(v)  to  zero.  We  shall  count  such  a 
push  against  the  most  recent  event  before  the  push  that  made  e(v)  nonzero.  We  call  such  an 
event  an  activation  of  v.  We  shall  derive  a  bound  on  the  number  of  activations  of 
0{nm  +  in^/l)  log  t/)plus  0{nlz)  per  wave,  thereby  proving  the  lemma. 

Initializing  the  preflow/at  the  beginning  of  the  entire  algorithm  causes  0{n)  vertex  activa- 
tions. Every  other  activation  is  due  to  a  tree  push  (case  2b  of  tree-push' relabel).  If  such  a  tree 
push  results  in  a  link  or  cut,  we  charge  the  activation  against  the  corresponding  link  or  cut.  Simi- 
larly, if  the  tree  push  is  saturating,  we  charge  the  activation  against  the  arc  saturation.  Such 
charges  account  for  0(nm)  activations. 

Any  remaining  activation,  say  of  a  vertex  x,  is  produced  by  a  tree  push  through  an  arc  (v,h) 
to  X  =  find-root{w),  after  which  find-size(v)  +  find-sizeiw)  >  z,  since  no  link  is  performed.  Let  T^ 
and  T^  be  the  dynamic  trees  containing  v  and  w,  respectively.  We  call  a  dynamic  tree  large  if  it 
contains  more  than  z/2  vertices  and  small  otherwise.  One  of  Ty  and  r»  must  be  large  just  after 
the  push. 
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If  the  tree  push  is  big,  we  charge  the  activauon  of  x  against  the  push.  By  the  argument 
above  counting  big  tree  pushes,  this  accounts  for  0(n^  +  (n^//)  log  U)  activations.  Otherwise, 
the  push  is  little,  and  thus  it  reduces  e(v)  to  zero.  There  can  be  at  most  one  such  push  from  v  per 
wave.  If  Ty  has  changed  between  the  beginning  of  the  most  recent  wave  and  the  push,  we  charge 
the  activation  of  x  to  the  link  or  cut  that  most  recently  changed  T^.  There  are  Oinm)  such  charges, 
at  most  one  per  link,  at  most  two  per  cut,  since  vertex  v  is  the  root  of  Ty  when  the  push  occurs. 
Similarly,  if  T^  has  changed  between  the  beginning  of  the  most  recent  wave  and  the  push,  we 
charge  the  activation  of  x  to  the  link  or  cut  that  most  recently  changed  T^.  There  are  0(nm)  such 
charges,  since  a  vertex  x  can  be  activated  at  most  once  per  wave,  and  x  is  the  root  of  T^  when  the 
push  occurs.  If  neither  Ty  nor  r»  has  changed,  we  charge  the  activation  of  x  to  whichever  of  Ty 
or  T^  is  large.  Each  large  tree  existing  at  the  beginning  of  a  wave  can  be  charged  at  most  twice 
(once  as  a  Ty,  once  as  a  r»).  Since  tfiere  are  fewer  than  2n/z  large  trees  at  the  beginning  of  any 
wave,  the  total  number  of  such  charges  is  0(n/z)  per  wave.  Combining  all  our  estimates  gives 
the  claimed  bound  on  vertex  activations.  D 

Theorem  5.4.  The  number  of  tree  pushes  done  by  the  tree  scaling  algorithm  is  0(nm)  if  z  is 

chosen  equal  to  min  (n,  max  (1, — :r  logL'})  and  /  is  chosen  equal  to  m/n  if  z  =  1,  —  logU  if 

m^  m 

\  <z<n. 

Proof.  By  Lemmas  5.2  and  5.3,  the  number  of  tree  pushes  done  by  the  algorithm  is 
Oinm  +  in^/l)  logU  +  —  min  [nk  +  logt/.n^ )).  We  consider  three  cases; 

Case  1:      — 7-logf/Sl.       Then  z  =  l   and  l  =  mln.    The  number  of  tree  pushes  is 
m 

0{run  +  {n^ /m)\ogU  +  n  minim  +  ]ogU,  n^]).  Since  logU  ^m^/n^,  the  bound  on  tree  pushes 
is  Oinm -i-  min{m^/n,  n^})  =  0(nm),  because  Tn}ln<,n}  implies  m<ji^,  which  means 
rn^ln  ^  nm,  and  m}ln  S  n}  implies  m  t  n^,  which  means  n'  ^  nm.  (This  analysis  allows  for  the 
possibility  of  multiple  arcs  in  G.) 


„2  f^l                            „ 

Case  2:     1  <  — r-logU  <  n.  Then    z  =  — ylogU    and    /  =  — \ogU.     The  number  of 

2  2                                                                                                 2 

tree  pushes  is    C>(n/rj+( — ; )  min{ — logl/ +  logt/,n^))  =  0(nm  +  min  (m^/n.- — —])  = 

n  logty  m                                                                  iogu 

0{nm  +  min{m^/n,  n^)),  since  \ogU  <m^/n^.   Since  min(m^/n,fi')  ^  nm  (see  Case  1).  the 
number  of  tree  pushes  is  Oinm). 


17 


^2  „ 

Case  3:     — r-  \ogU  ^  n.     Then  z  =  n  and  /  =  —  logU.      The  number  of  tree  pushes  is 


0(wn  +  min( — logU  +  logU.n^ ])  =  0(.nm).  D 
m 


f^2 
Observe  that  if  z  is  chosen  as  in  Theorem  5.4,  i.e.  equal  to  min{n,  max  {1, — r-  \ogU] ),  the 

m 

time  per  dynamic  tree  operation  is  0(\og(z  +  1))  =  0(\og  ( —  Qogf/)"^)).  Thus  the  running  time 

m 

of    the     tree     scaling     algorithm,     with     k     and     z    chosen     as     in    Theorem     5.4,     is 

0(nm  log  ( —  Gogt/)"^  +  2))  (igrwring  time  spent  manipulating  L). 
tn 

All  that  remains  is  to  show  that  the  manipulations  of  the  list  L  in  wave  steps  can  be  per- 

fonmed  in  0{nm  log( — Gogt/)"^  +  2))  time.  It  is  not  sufficient  to  represent  L  simply  as  a  lirJced 
m 

list,  for  then  the  time  scanning  L  will  be  0{n)  per  wave,  and  this  scanning  time  will  dominate  the 
time  for  the  rest  of  the  computation.  Instead,  we  represent  L  as  a  list  of  sublists,  each  of  which  is 
in  turn  represented  by  a  balanced  tree.  This  idea  was  used  by  Goldberg  and  Tarjan  in  a 
minimum-cost  flow  algorithm  [11].  Our  version  of  the  method  differs  from  theirs  in  that  we 
represent  the  sublists  by  ordinary  balanced  search  trees  instead  of  finger  search  trees,  and  we  limit 
the  size  of  each  sublist  to  z.  the  maximum  tree  size  in  the  tree  scaling  algorithm.  These  changes 
lead  to  an  analysis  that  parallels  the  analysis  of  the  tree  pushes. 

The  details  of  the  data  structure  are  as  follows.  The  list  L  is  maintained  throughout  the 
computation  so  that  its  vertices  are  in  topological  order  with  respea  to  the  eligible  arcs;  this  elim- 
inates the  need  to  reconstruct  L  at  the  beginning  of  each  wave.  Maintaining  L  in  tofwlogical 
order  merely  requires  moving  each  relabeled  vertex  to  the  front  of  L  when  it  is  relabeled. 

List  L  is  maintained  as  a  doubly-linked  list  of  sublists.  These  sublists  have  two  properties: 


(a)  all  sublists  have  size  at  most  z; 

(b)  any  active  vertex  is  last  on  its  sublisL 

The  sublists  are  represented  as  balanced  search  trees,  e.g.  red-black  trees  [20]  or  splay  trees 
[19,20].  This  representation  supports  the  follosving  op>erations  on  the  sublists: 

(i)  Initialize  a  sublist  to  contain  a  given  single  vertex. 
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00         Find  the  number  of  vertices  in  a  given  sublisl. 

Oii)        Find  the  sublist  containing  a  given  vertex. 

Ov)        Find  the  last  vertex  in  a  given  sublist. 

(v)         Split  the  sublist  containing  a  given  vertex  just  before  (or  just  after)  the  vertex  and  return 
the  two  resulting  sublists. 

(vi)        Concatenate  two  sublists  and  return  the  resulting  sublist. 

The  time  required  for  these  operations  is  0(1)  for  each  operation  of  the  types  (i)  and  (ii), 
and  0(\og(z  +  1))  for  each  operation  of  types  (iii)-(vi),  assuming  that  no  sublist  size  ever  exceeds 

2. 

During  the  tree  scaling  algorithm,  the  computations  done  on  L  are  as  follows.  Initially, 
each  vertex  other  than  s  and  /  is  initialized  to  be  a  single-vertex  sublisl,  and  these  sublists  are 
ordered  in  nonincreasing  order  by  vertex  label.  Each  time  a  vertex  v  is  activated  (i.e.  e(v) 
becomes  nonzero),  the  sublist  containing  v  is  split  just  after  v.  This  maintains  property  (b). 
When  a  vertex  v  is  relabeled,  the  sublisl  containing  v  is  split  just  before  v,  and  the  (singleton)  sub- 
lisl containing  v  is  moved  from  its  current  position  to  the  front  of  L. 

Concatenation  of  sublists  takes  place  only  during  waves.  A  wave  consists  of  initializing 
vertex  v  to  be  the  last  vertex  on  the  first  sublist  of  L,  and  repeating  the  following  step  until  v  is 
null: 

Scan.  Apply  the  appropriate  one  of  the  following  cases: 

Case  1 :  V  is  active.  Perform  stack-push/relabel{v),  splitting  sublists  at  newly  active  ver- 
tices and  moving  relabeled  vertices  (except  for  v)  to  the  front  of  L  as  described  above.  If  v 
has  been  relabeled,  let  j:  be  v,  and  proceed  as  follows:  Replace  v  by  the  last  vertex  on  the 
sublist  after  v,  or  by  null  if  there  is  no  such  sublisl;  split  x  (the  old  v)  from  its  sublist  and 
move  j:  to  the  front  of  L. 

Case  2:  v  is  inactive.  If  v  is  in  the  last  sublisl,  replace  v  by  null.  Otherwise,  let  B  and  B' 
be  the  sublists  contairung  v  and  the  vertex  after  v,  respectively.  Replace  v  by  the  last  ver- 
tex in  fl'.  If  B  and  B'  together  contain  no  more  than  z  vertices,  concatenate  B  and  B'. 
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Note  that,  during  a  wave,  if  a  stack-pushJ relabel  applied  to  a  vertex  v  does  not  relabel  v, 
then  V  is  immediately  reprocessed  in  Case  2  oi  scan. 

The  analysis  of  the  list  processing  that  takes  place  during  waves  is  much  like  the  analysis 
already  done  to  count  tree  pushes.  The  time  spent  processing  L  during  a  execution  of  case  1  is 
O  Gog  (z  +  1))  per  tree  push  plus  O  Gog  (z  +  1))  per  relabeling.  The  time  spent  processing  L  dur- 
ing an  execution  of  case  2  is  O  Gog  (z  +  1)).  Except  for  O  (n/z)  executions  of  case  2  per  wave 
(those  in  which  B  and  B '  together  contain  more  than  z  vertices),  each  execution  of  case  2  causes  a 
sublist  concatenation.  The  number  of  concatenations  can  exceed  the  number  of  splits  by  at  most 
n-1,  and  the  number  of  splits  is  at  most  the  number  of  tree  pushes  plus  the  number  of  relabelings. 
Thus  we  obtain  the  following  result: 

Lemma  5.5.    The    total   time    spent   processing   L   during   the    tree    scaling    algorithm    is 
0{nm  log  (z  +  1)),  if  z  and  /  are  chosen  as  in  Theorem  5.4. 

Proof.  The  number  of  relabelings  is  O  (nm).  So  is  the  number  of  tree  pushes,  by  Theorem  5.4.  D 


Theorem  5.6.   The  tree  scaling  algorithm  runs  in  O  {nm  log  ( —  Gog  U)      +2))  time  if  z  is 

m 

2 

chosen  equal  to  min{rt,  max{l, — r-  log  U}  and  /is  chosen  equal  to  —  if  z  =  1  or  to  —  log  f/ if 

m^  n  m 

2>  1. 


Proof.         Immediate        from        Theorem        5.4        and        Lemma        5.5,        since 
log  (z  +  1)  =  O  Gog  (—  Gog  Uf^  +  2)). 


Remark.  The  bound  in  Theorem  5.6  can  be  improved  to  O  (nm  log  (        ^ —  +  2),  analo- 


\I/2 

m 

gously  to  the  improvements  in  the  bounds  of  Sections  3  and  4.  To  obtain  this  improvement,  we 
run  the  stack  scaling  algorithm  until  A  ^  U*,  and  then  switch  to  the  tree  scaling  algorithm.  D 


6.  Bounds  in  a  Semilogarithmic  Computation  Model 

The  time  bounds  derived  in  Sections  3-5  are  all  based  on  the  assumption  that  addition  and 
comparison  of  integers  of  magnitude  at  most  U  takes  0(1)  time.  If  U  is  huge,  this  assumption 
may  be  unjustified.  In  this  section  we  explore  the  consequences  to  our  results  of  using  a  semilog- 
arithmic computation  model,  in  which  each  computer  word  is  allowed  to  hold  O  (logn)  bits  and 
any  operation  involving  0(1)  words  takes  0(1)  time.  In  this  model,  all  the  elementary  graph  and 
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list  manipulation  operations  needed  by  our  algorithms  take  0(1)  time  each,  but  adding  or  com- 
paring two  capacity  or  flow  values  can  take  O  Oog„  V)  time  if  an  exact  answer  is  required.  Thus 
a  straightforward  translation  of  our  results  into  the  semilogarithmic  model  increases  each  of  the 
time  bounds  by  a  factor  of  0  Gog„  t/). 

By  suitably  modifying  the  algorithms,  however,  it  is  possible  to  reduce  the  extra  time  each 
of  our  algorithms  requires  in  the  semilogrithmic  model  to  an  additive  term  of  O  (m  log,  U). 
Thus     for     example     the     running     time     of     the     tree     scahng     algorithm     becomes 

O  (nm  log  ( —  Gog  f/)"^  +  2)  +  m  log,  U).  Note  that  the  time  needed  to  read  all  the  arc  capaci- 
m 

ties  is  6  (m  log,  U)  in  this  model. 

The  general  approach  is  to  approximately  solve  a  sequence  of  O  Qog,  U)  closer-and-closer 
approximations  to  the  original  problem.  Each  approximation  uses  as  a  starting  point  the  approxi- 
mate solution  to  the  previous  problem.  Solving  each  problem  requires  manipulation  of  integers 
of  only  O  Gog  n)  bits. 

Making  this  approach  work  involves  a  number  of  messy  technical  detaQs.  Since  the  result  is 
mainly  of  theoretical  interest,  we  shall  merely  sketch  the  ideas  involved  in  modifying  the  algo- 
rithms of  Sections  4  and  5.  We  assume  (without  loss  of  generality)  that  n  is  a  power  of  two  and 

that  t/  is  a  power  of  n.  For  k  =  1,2 log,  U,  we  define  problem  )t  to  be  the  maximum  flow 

problem  on  the  graph  G  with  arc  capacities  w*  (v,^)  defined  by  Ut(v,w)=\  u(v.u')/5tl  5^, 
where  5i  =  t///i*.  Observe  that  for  k  =  log,  t/,  5*  =  1;  thus  problem  log,  f/  is  the  original  prob- 
lem. Also,  for  a  general  value  of  k,  all  capacities  in  problem  k  are  divisible  by  5*;  and.  for  every 
arc  (v.w),  0  ^  «t_i  (v.w)  -Uk(v,w)  ^  6i_i . 

We  solve  problem  kfoTk  =  1.2, ....log,  (7-1  approximately  and  then  solve  problem  log,  U 
exactly.  Throughout  the  computation  we  maintain  a  preflow  /  and  a  valid  labeling  d  for  the 
current  problem.  Preflow /satisfies  the  following  constraint: 


e{v)  ^  min  {     "^     f(u,v),  n(5t-\) )  for  all  v  *^  (flow  holdback  constraint).  (5) 

/(y.v)  >  0 


(In  inequality  (5).  the  sum  over/(u,v)  is  taken  to  be  zero  if /(u,v)  S  0  for  all  arcs  (u,v).) 

For  a  vertex  v,  the  available  excess  at  v.  which  is  the  amount  actually  allowed  to  be  pushed 
from  V,  is  a(v)  =  max  {0,e(v)-n(5i-l) ).  Throughout  the  computation,  the  available  excesses  are 
used  in  place  of  the  actual  excesses  to  determine  when  vertices  are  active  and  how  much  flow  can 
be  pushed. 

The  flow  holdbacks  are  needed  to  maintain  the  preflow  property  when  converting  a  solution 
to  one  problem  into  a  good  initial  preflow  for  the  next  problem.  For  the  first  problem,  the  preflow 
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/and  valid  labeling  d  are  initialized  exactly  as  in  Section  2  (using  the  arc  capacities  «i(v,w)). 
Once  a  solution  to  problem  k-l  is  computed,  it  is  converted  into  an  initial  preflow  for  problem  k 
by  replacing  the  current  preflow /by  the  preflow/'  defined  by /'(v,H')  =  min{/(v,H'),M4(v,>v)). 
Since  0^  M4_i(v,H')-«t((v,H')  <  5i_i,/'  is  obtained  from  /by  decreasing  the  flow  on  each  arc 
(v.w)  by  an  amount  between  zero  and  mm[f{v,w),  5t-i )  (inclusive).  Since  S^.j  =  nS*.  the  vali- 
dity of  the  flow  holdback  constraint  for /in  problem  k-l  implies  the  validity  of  the  flow  holdback 
constraint  forf  in  problem  k.  The  flow  holdback  constraint  implies  that  f  is  a  preflow,  since  f 
satisfies  the  capacity  constraint  by  construction.  Since  every  arc  saturated  by  /  in  problem  k-l 
remains  saturated  by  f  in  problem  k.d  is  a  valid  labeling  for/  in  problem  k. 

For  I  <k  <  log,  U,  a  preflow  / and  valid  labeling  d  constitute  a  solution  to  problem  k  if 
a(v)  <  2n(5i-l)  for  every  active  vertex  v.  Thus  a  solution  to  problem  log,  U  gives  a  maximum 
flow  in  the  original  network.  After  the  initialization  of  the  preflow  for  problem  k,  every  active 
vertex  has  an  excess  of  at  most  4n^5t.  of  which  3n^5t  =  3n5t_i  comes  from  the  excess  on  v  in 
problem  k-l  (since  a(v)  <  2n(5t_i-l)),  and  n^5t  =  n5k-i  comes  from  the  changes  made  to /in 
initializing  the  preflow  for  problem  k.  (These  are  overestimates.)  The  initial  value  of  A  (the 
bound  on  maximum  available  excess)  for  problem  k  is  4n^5k.  Problem  k  is  solved  by  performing 
a  number  of  phases,  continuing  until  the  maximum  available  excess  is  at  most  2n  (5t-l),  which 
happens  by  the  time  A  is  reduced  to  nS*  (or  to  1/2  if  /t  =  log,  U).  U  k  <  log,  U,  the  number  of 
phases  needed  is  at  most  log  n  +  2;  if  k=  log,  U,  the  number  of  phases  needed  is  at  most 
21ogn  -I-  3.  Thus  each  problem  requires  OQogn)  phases,  and  the  total  number  of  phases  over  all 
subproblems  is  O  (logn  log,  U)  =  0  (log  U). 

During  the  solution  of  problem  k,  the  preflow /is  maintained  so  lhat/(v,H')  is  a  multiple  of 
5i  for  each  arc  (v.w).  (This  happens  automatically,  since  A.  the  initial  flow  values,  and  all  capa- 
cities are  multiples  of  5^).  Thus  the  flow  values  and  capacities  can  be  represented  in  units  of  6*. 
Furthermore,  the  flow  values  and  capacities  are  not  represented  explictly,  but  rather  as  differences 
from  the  initial  flow  values.  To  be  more  precise:  let/4_i  be  the  final  preflow  for  problem  k-l\ 
let  /o  =  0.  Then  the  current  preflow  /is  represented  by  the  value  f{v,'w)-fk-\  (v.w)  for  each  arc 
(v,^).  The  capacity  fuiKtion  u^  is  represented  by  the  value  min  {«t(v,H')-/t_i(v,H'),  9n'  S^)  for 
each  arc  (v.w).  We  claim  that  the  amount  by  which  the  flow  on  an  arc  can  change  during  the  solu- 
tion of  problem  k  and  subsequent  problems  is  at  most  8/j'  5*.  This  claim  implies  that  if 
"i  (v.H')-/t(v,H')  >  9n'  5*,  the  extra  residual  capacity  (beyond  9n'  5*)  on  arc  (v.w)  will  never  be 
used  and  can  be  ignored.  To  prove  the  claim,  we  note  that  between  relabelings  at  most 
n  A  <4n'  5*  units  of  flow  can  be  moved  through  any  given  arrc,  and  there  are  at  most  In}  rela- 
belings over  the  entire  algorithm.  This  accounts  for  8n'  5*  units  of  change.  The  change  on  an 

arc  due  to  modifying  /between  subproblems  is  at  most    J^    6^  <  2hk.\  =  2n  5*.  Thus  the  claim 

j-k-\ 
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is  true. 

Storing  /and  ut  in  such  a  difference  form,  in  units  of  6*,  allows  the  algorithm  to  manipulate 
numbers  consisting  of  only  0(\ogn)  bits,  and  all  necessary  additions  and  comparisons  of  flows 
and  capacities  take  0(1)  time.  Constructing  the  final  preflow  is  merely  a  matter  of  adding 
together  the  successive  flow  differences /i-/o,/2-/i.--  computed  when  solving  successive  prob- 
lems. By  adding  these  differences  from  right  to  left,  the  time  per  addition  can  be  made  0(m), 
and  the  total  time  is  0(m  log^  U).  This  is  also  the  total  time  required  for  initializing  the  preflow 
for  each  problem  (Oim)  time  per  problem). 

If  the  wave  scaling  algorithm  is  used,  the  bounds  in  Section  4  remain  valid  almost  without 
modification,  and  apply  to  the  calculations  involved  in  solving  all  log„  U  problems.  Specifically, 
the  number  of  relabelings  is  O(n^),  the  number  of  saturating  pushes  is  0(nm),  the  number  of 
nonsaturaiing  pushes  is  O  (n^  +  (n^/l)  log  U)  plus  0(n)  per  wave,  and  the  number  of  waves  is 
0(min  [nl  +  log  U,  n^  +  log,  U]).  Tht  bound  on  waves  is  the  only  one  that  changes,  in  that  n^ 
becomes  n^  +log;,  U;  this  is  because  a  wave  without  relabelings  terminates  only  the  solution  of 
the  current  problem,  and  rwt  the  entire  algorithm.  Chooseing  /  =  Gog  U)^'^  gives  an  overall 
bound  of  0{nm  +  n^  (log  l/)"^  +  m  log,  if)  time. 

Use    of   the    tree    scaling    algorithm    of   Section    5    results    in    a    time    bound    of 

Oinm  log( — Gog  t/)"^  +  2)  +  m  log,  U).  To  obuin  this  bound,  we  must  verify  that  the  extra 
m 

time  required  for  initializing  the  list  L  and  the  dynamic  trees  for  each  problem,  and  the  extra  time 
for  tree  pushes  associated  with  vertices  that  become  active  at  the  beginning  of  a  new  problem,  is 
0(m)  per  problem. 

We  need  two  extra  facts  about  the  dynamic  tree  data  structure.  By  traversing  the  entire 
data  structure,  the  flow  values  for  all  the  dynamic  tree  arcs  can  be  computed  in  0{n)  time.  Furth- 
ermore, a  new  set  of  flow  values  for  all  these  arcs  can  be  installed  in  the  data  structure  in  0(n) 
time.  These  facts  can  easily  be  verified  by  checking  the  description  of  the  data  structure  [19.20]. 

When  the  flow /is  modified  at  the  beginning  of  a  new  problem,  new  vertices  with  available 
excess  are  created,  possibly  as  many  as  n-2.  Handling  these  newly  active  vertices  is  the  hardest 
part  of  the  initialization  of  the  new  problem.  To  begin  a  new  problem,  the  current  flow  on  all 
dynamic  tree  arcs  is  computed  explicitly.  (This  takes  0(n)  time  as  noted  above.)  Then  the  flow 
and  capacities  are  modified  to  their  initial  values  for  the  new  problem.  The  bound  on  maximum 
available  excess  is  set  equal  to  4n'  6*.  (This  is  n  times  the  value  proposed  eariier  in  this  section, 
and  is  greater  than  the  total  excess  on  all  active  nodes.) 

Next,  the  vertices  of  C  are  scanned  in  topological  order  with  respect  to  the  set  of  eligible 
arcs.  A  vertex  v  is  scanned  by  applying  to  it  a  pushlrelabel  step  of  the  kind  defined  in  Section  2. 
Such  a  step  either  relabels  v  or  reduces  its  available  excess  to  zero.  After  all  vertices  have  been 


23 


scanned,  all  the  available  excess  is  on  vertices  that  have  been  relabeled  during  the  scanning.  The 
current  flow  values  for  dynamic  tree  arcs  are  reiosuied  in  the  d>-namic  tree  data  structure,  and 
every  dyiiamic  tree  arc  that  has  been  saturated  during  the  scanning  or  that  has  had  one  of  its  end 
vertices  relabeled  is  cut  Finally,  the  list  L  is  reinitialized  by  sorting  the  vertices  in  topological 
order  with  respect  to  the  set  of  eligible  arcs  and  constructing  one  sublisi  in  L  for  each  active  ver- 
tex V.  with  V  as  its  last  vertex,  and  one  sublist  for  any  vertices  in  L  that  foUov^'  the  last  active  ver- 
tex. Constructing  all  the  sublists  takes  0{n)  time,  since  a  binar>-  search  tree  containing  x  items 
can  be  initialized  in  0{x)  time  by  doing  repeated  pair-wise  concatenanon  (concatenating  pairs  of 
singletons,  then  pairs  of  pairs,  and  so  on).  Then  the  tree  scaling  algorithm  is  begun. 

Because  the  starting  value  of  A  for  each  problem  is  increased,  the  number  of  phases 
increases,  but  only  by  \ogn  per  subproblem,  or  logl^  overall.  The  time  to  carry  out  all  the  scan- 
ning in  the  initialization  is  0{m).  The  time  to  initialize  the  d>Tiamic  tree  data  structure  is  0{jn) 
plus  (9Gog(r  -K  1))  per  cut;  the  time  for  a  cut  can  be  charged  against  the  corresponding  arc  satura- 
tion or  vertex  relabeling.  The  time  to  initialize  L  is  0{m).  Once  the  initialization  has  been  com- 
pleted, the  only  active  vertices  are  those  that  have  been  relabeled  during  the  initialization,  and  the 
timing  analysis  for  the  tree  scaling  algorithm  is  virtually  the  same  as  the  analysis  in  Section  5. 

Thus  we  obtain  a  total  time  bound  of  O  (nm  log  ( — Gog  t/)"^)  +  2)  +  m  log^,  U). 

m 

The  techniques  discussed  above  can  also  be  applied  to  the  maximum  flow  algorithm  of 

Goldberg  and  Tarjan   [8,10],  resulting  in  a  time  bound  in  the  semilogarithmic  model  of 

O  {run  log(/i^/m)  -fm  log^,  U).  The  details  are  straightforward. 


7.  Remarks 

The  algorithms  of  Section  3  and  4  not  only  have  a  good  theoretical  time  bounds,  but  they 
may  be  very  efficient  in  practice.  We  hope  to  conduct  experiments  to  determine  whether  either  of 
these  algorithms  is  competitive  with  previously  known  methods.  The  tree  scaling  algorithm  of 
Section  5  is  perhaps  mainly  of  theoretical  interest,  although  for  huge  nenAorks  there  is  some 
chance  that  the  use  of  dvnamic  trees  may  be  practical. 

The  two  obvious  open  theoretical  questions  are  whether  further  improvement  in  the  time 
bound  for  the  maximum  flow  problem  is  possible  and  whether  our  ideas  extend  to  the  minimum 
cost  flow  problem.  It  is  not  unreasonable  to  hope  for  a  bound  of  0{nm)  for  the  maximum  flow 
problem;  note  that  the  bounds  of  Theorems  4.3  and  5.6  are  0{.nm)  for  graphs  that  are  not  too 
sparse  and  whose  arc  capacities  are  not  too  large,  i.e.  Oogf.')'  *  =  0{m  'n).  Obtaining  a  bound 
better  than  0{nm)  would  seem  to  require  major  new  ideas. 
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As  a  pania]  answer  to  the  second  question,  we  have  been  able  to  obtain  a  time  bound  of 
0(nm  loglog  U  log  (nC))  for  the  minimum  cost  flow  problem,  where  C  is  the  maximum  arc  cost, 
assuming  all  arc  costs  are  integral  [2].  This  result  uses  some  ideas  in  the  present  paper  and  some 
additional  ones. 
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